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Abstract. We use the deepest and most complete redshift catalog currently avail- 
able (the Hubble Deep Field (HDF) North supplemented by new HDF South 
redshift data) to minimize residuals between photometric and spectroscopic red- 
shift estimates. The good agreement at Zspcc < 1.5 shows that model libraries 
provide a good description of the galaxy population. At 2;spcc > 2.0, the system- 
atic shift between photometric and spectroscopic redshifts decreases when the 
modeling of the absorption by the interstellar and intergalactic media is refined. 
As a result, in the entire redshift range z G [0,6], residuals between photometric 
and spectroscopic redshifts are roughly halved. For objects fainter than the spec- 
troscopic limit, the main source of uncertainty in photometric redshifts is related 
to photometric errors, and can be assessed with Monte Carlo simulations. 
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1. Introduction 

Reliable spectroscopy of galaxies fainter than /ab ~ 24 is still an unattainable goal even 
with new-generation telescopes. To estimate the redshift for a majority of these objects 
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one must therefore rely on multicolor observations in order to reconstruct, at very low 
resolution, the apparent spectral energy distribution (SED). 

Much efforts have been devoted, in the recent literature, to improve the accuracy 
of the photometric approach compared with direct spectroscopy. Photometric redshift 
techniques have been extensively applied, in particular, to the ultra-deep images of the 
Hubble Deep Field North (HDFN, Williams et al. 1996, see also Ferguson, Dickinson and 
Williams 2000) to study the luminosity function, the evolution of the star formation rate 
(SFR) density, and the clustering properties of the galaxies observed there (Madau et al. 
1996; Connolly et al. 1997; Sawicki et al. 1997; Franceschini et al. 1998; Pascarelle et al. 
1998; Connolly et al. 1998; Miralles & PeUo 1998; Roukema et al. 1999, Arnouts et al. 
1999). 

The sample of relatively bright HDFN objects with reliable spectroscopic redshifts is 
steadily increasing and now amounts to about 150 objects, a value that nearly doubles the 
initial list of Sawicki et al. (1997). About 20% of these objects are beyond z > 2.0, thus 
allowing a more robust statistical comparison between spectroscopic and photometric 
estimates of z even for the most distant objects in the field. 

We previously discussed (Massarotti et al. 2001a, Paper I hereafter) the performance 
of the SED fitting method as a function of the template set adopted to reproduce galaxy 
colors. Our results suggested that, as long as both normal and starburst galaxies are 
taken into account in the reference templates, the photometric redshift at z < 1.5 is 
constrained by a reliable detection of the 4000 A Balmer break in the target SED. At 
larger redshifts, on the contrary, the main constraint is set by modeling the effects of 
the interstellar (ISM) and intergalactic (IGM) absorption, which modulate the SEDs at 
different wavelengths. In this paper, we further discuss some important aspects of this 
problem, and suggest an improved model for the ISM and IGM absorption which sensibly 
increases the accuracy of the method at z '^2. 

Photometric redshifts are estimated by comparing the observed broad-band colors of 
galaxies with the library of model fluxes provided by the code of Buzzoni (1989, 1995, 
2001; hereafter BUZ). Results obtained with the codes of Bruzual & Chariot (1993, BC) 
and Fioc & Rocca-Volmerange (1997, FRV) are also reported, since these theoretical 
models have also been extensively used in the literature. 

The key steps of our procedure are summarized in Sect. 2, where we also introduce 
the model libraries. The redshift database used in our analysis is presented in Sect. 3, 
comparing the ^;spec distribution with our ^;phot output. The ISM and IGM effects are 
more specifically dealt with in Sect. 4, while in Sect. 5 we discuss the match between 
models and observations even at the faintest magnitude limit of the HDFN data. Section 
6 summarizes our conclusions. 
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2. Input ingredients and fitting algorithm 

As in Paper I, the photometric catalog of Fernandez-Soto et al. (1999) will be used as a 
reference in our analysis. It is made up by 1067 objects identified from the HST WFPC2 
images in the four photometric bands, namely f/aoo, -B450, Veoe, and /814, and completed 
with infrared observations in the Johnson J, H, and K bands with the IRIM camera at 
KPNO (Dickinson 1998). 

Galaxy model libraries have been assembled following the strategy of Paper I. A 
notable change in the present calculations involves a slight enlargement of the set of 
starburst models in order to cover a wider range of physical parameters and ensure a 
more accurate description of the colors of the most distant galaxies. 

To the Leitherer et al. (1999) starburst sequence of Paper I, namely for t = 50, 100, 
500, and 800 Myr, with Zq and Salpeter IMF, we added here two supplementary models 
at f = 10 and 250 Myr, and included evolution for Zq/20 metallicity in addition to the 
solar case previously considered. 

In order to account for recent star formation superimposed on quiescent red (elliptical) 
galaxies, we also considered a new set of templates. These are obtained adding the 50 
Myr starburst spectrum (with Salpeter IMF and Zq) to simple stellar population models 
for a multi-component fit to the galaxy SED. In this way the overall color distribution 
of early type galaxies at ^; < 1.5, as found by Franceschini et al. (1998) and Rodighiero 
et al. (2001), is much better accounted for. 

As a first step, the ISM and IGM absorptions have been modeled in as in Paper I. 
Galaxy model fluxes have been corrected for ISM effects following Calzetti (1999). As in 
Paper I, the E{B — V) color excess has been set as a free fitting parameter which could 
take the discrete values 0.0, 0.05, 0.1, 0.2, 0.3, 0.4. The IGM absorption as a function of 
redshift follows the analaysis by Madau (1995). 

According to our SED fitting method, we searched for the "best template" solution as 
which minimizes the general function (see Eq. (1) of Paper I). In the multi-component 
fit to the galaxy SED the functional takes the form: 



burst model, respectively, scaled by Sgsp and Sburst, and is the number of filters. 
For galaxies with Zgpec available (A'g) we define two simple statistical moments: 
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to study residuals between photometric and spectroscopic redshift values. 

Catastrophic outliers (that is, galaxies with l^spcc ~ -ZphotI > 1-0) are clipped from the 
galaxy distribution to compute representative values for Az and ct^- 

3. The spectroscopic sample 

Cohen et al. (2000) provided an exhaustive compilation of all spectroscopic redshifts 
measured with the Keck telescope in the HDFN proper and flanking fields. The work 
includes previously published data, new redshifts measurements, and corrections to erro- 
neous values given in the literature. This list represents the deepest and most complete 
redshift catalog for the HDFN to date. 

Fernandez-Soto et al. (2001) matched photometric and redshift catalogs in the HDFN 
proper. They cross-correlated coordinates and magnitudes, resulting in a list of 146 
sources (see Table 2 therein). We adopt here their naming convention and follow their 
considerations regarding some galaxies with discordant Zgpcc vs. Zphot values, except for 
the following few objects. 

HDr36441— 1410: we assume that the spectroscopic measure Zgpcc = 2.267 is cor- 
rect. Fernandez-Soto et al. (2001) claimed that features identifications in the spectrum 
are compatible with their estimated 2;phot = 0.01, but HDF36441— 1410 is a [/ drop- 
out star-forming galaxy ([/300 — S450 ~ 1.7, Vgoe — -^814 ~ 0.0). These colors (blue in 
^606 ~ -^814, red in Usoa — S450) can be explained only as consequence of the IGM ab- 
sorption of UV photons in the C/300 band, and are incompatible with z ^ starburst 
colors. 

HDF36446— 1227: we prefer the spectroscopic value Zgpcc — 2.035 proposed by 
Lowenthal et al. (1997), to the (revised) value Zgpcc = 2.500 reported by Cohen et al. 
(2000) and accepted by Fernandez-Soto et al. (2001). In the redshift range Zgpcc G [2, 2.3] 
there is no galaxy bluer than HDF36446— 1227 in C/300 ~ -S450J and according to the 
starburst templates of our models its optical colors are incompatible with Zgpoc > 2.3 
(see Sect. |4^ ). 

HDF36396 1230, 36494 1317, and 36561 1330: for these three galaxies we fol- 
low Cohen (2001), who adopts the redshift values given by Fernandez-Soto et al. (2001). 

HDF36478— 1256: with Zgpcc — 2.931, this galaxy was tagged as a catastrophic 
outlier in all previous works (see e.g. Fernandez-Soto et al. 1999; Bolzonella et al. 2000; 
Fernandez-Soto et al. 2001; Fontana et al. 2000). Arnouts et al. (1999) suggest that 
the photometry could be incorrect because of the complex source morphology. Given 
the unclear situation and the large discrepancy for this galaxy (our photometric redshift 
estimate is Zphot — 0.26) we excluded it from the catalog. 

Six new galaxies from the HDF-South (HDFS) have also been added to our list. 
Three of them are from Glazebrook et al. (2001) with Zgpoc — 0.580, 0.582, and 0.696, 
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Fig. 1. Comparison between spectroscopic and photometric redshifts according to differ- 
ent template libraries (see the text for acronyms). The solid line is for Az = 0.0. Note 
the different scale for panel d. 



and three are from Cristiani et al. (2000) with Zspcc = 0.565, 2.79, and 3.2 respec- 
tively. The photometry was obtained from the web page of the Stony Brook group (see 



http:/ /www. ess. sunysb.edu/astro/hdfs ) in the four WFPC2 filters, and in the J, i?, and 



K infrared bands by the SOFI camera at NTT (Da Costa et al. 2001). 

Our final spectroscopic catalog contains 146 objects, 34 of them with Zspoc > 2. 



3.1. Comparison between z^pec snd Zphot 

In Fig. |l| we compare spectroscopic and photometric redshifts obtained with different 
sets of reference galaxy models. 

At Zspoc < 1.5 (panels a,h,c) the agreement is excellent, regardless of the template 
set chosen. Dispersion values of the photometric redshift estimates are — 0.078, 0.074, 
0.072 for the BUZ, FRV and BC libraries, respectively, while systematic errors are Az = 
0.007, 0.026, 0.004. There are no galaxies with jzspcc — ^photl > 0.5. 

Since in this redshift range the accuracy of the method depends mainly on the ability 
of models to reproduce the 4000 A Balmer break along the whole Hubble morphological 
sequence (including starburst galaxies), one can conclude that all model libraries contain 
enough information to match properly the SEDs. 
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For the same sample, at Zspcc < 1-5, Fontana and collaborators obtained CTz = 0.084 



and Az = 0.033 (see the web page http://www.mporzio.astro.it/HIGHZ/HDF.html). 

Fernandez-Soto et al. (1999) obtained CTz ~ 0.116 and Az ~ —0.007, further improved 
to (Tz = 0.107 and Az = —0.004 after including starburst galaxies in their template li- 
brary (Fernandez-Soto et al. 2001). Bolzonella et al. (2000) obtained tiz ~ 0.082, using 
a smaller spectroscopic sample, that did not account for the new spectroscopic redshifts 
from Cohen and collaborators. 

The situation becomes more complex at Zgpcc > 2 (see panel d of Fig. The dis- 
persion increases to <7z — 0.177 and there are two objects with jzspcc — ■ZphotI > 0.5. In 
this redshift range almost all galaxies are best fit by starburst models, and we therefore 
focussed on the Leitherer et al. (1999) templates. 

Part of the larger scatter can certainly result from the fainter magnitudes (lower 
signal-to-noise ratio), but a systematic shift (Az = 0.156) is significant in panel d, and 
cannot be explained in terms of photometric errors only. The tendency to underestimate 
redshifts for Zgpoc > 2 seems to be a common feature of all photometric studies found in 
the literature, regardless of the specific choice of theoretical and/or empirical template 
library. 

The only exception in this regard are the results of Fernandez-Soto et al. (2001) 
which show systematics in the opposite sense, with Az ~ —0.136. These authors did 
not however take into account explicitely the ISM reddening modulation of the template 
SEDs (the average amount of dust that the empirical templates do contain may not be 
entirely representative of the dust content of galaxies at high redshifts). As we noted in 
Paper I, this causes a significant overestimate (by about 10-20%) of z at large distances, 
because the IGM absorption, which increases with redshift, is forced to mimic color 
reddening. 

The presence of dust in starburst galaxies at high redshift has been confirmed by 
analyzing galaxy colors with the SED fitting method (see Massarotti et al. 2001b, and 
references therein). Meurer et al. (1999) adopted an entirely empirical approach based 
on the correlation between color reddening, dust absorption and far infrared flux for 
local starburst galaxies (supposed to be also valid at high redshifts). According to their 
results (see also Meurer et al. 1997; Sawicki & Yee 1998), at < z >= 2.75 a color- 
luminosity correlation is already in place: more luminous galaxies (i.e. galaxies with 
available spectroscopic redshifts) have a larger dust content. 

The observed systematic shift in the Zspcc vs. Zphot distribution of Fig. |l|, calls therefore 
for a revised analysis of the ISM and IGM effects in the z 2 data. 
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4. A revised approach to IGM and ISM absorption 

In Paper I we showed that at high redshifts the SED fitting method provides stable 
estimates of galaxy redshifts, even when using an extremely poor template library, while 
the critical ingredients are dust reddening and IGM attenuation. It is important to explore 
in some detail the physical processes of photon absorption in both cases. 

4.1. The IGM opacity as a function of redshift 

Following Madau (1995), the Lya forest contribution to the attenuation of the SED 
caused by the IGM absorption can be written in terms of an effective optical depth as 

^intr(A) 

with 



(4) 



TcS 



A X (Aobs/Aa)i+^ (5) 



where Fobs (A, z) and Finti (A) are the observed and intrinsic fluxes respectively, Aobs = 
A X (1 + z), and A^ 1216 A. 

Press et al. (1993), in their analysis of the distribution of Lya forest lines in a sample 
of 29 quasars at z e [3.1, 4.5], found {A, 7) = (0.0036, 2.46), which are the values usually 
adopted in the literature. 

In a recent paper, Scott et al. (2000) studied a sample containing 97 quasars at 
z e [1.7,3.6], a redshift range complementary to that studied by Press et al. (1993). 
Considering absorption lines above the threshold W > 0.16 A, they obtain (A, 7) ~ 
(0.00759, 1.35). The IGM transmission at z < 3 is therefore larger than the one expected 
by extrapolating the Press et al. (1993) results at these redshifts. 

We decided to adopt the Scott et al. (2000) results to describe the IGM absorption 
at z < 3 and the Press et al. (1993) results at z > 4. In the intermediate redshift range 
(i.e. between 3 < z < 4), the A and 7 coefficients of Eq. (5) have been interpolated 
hnearly imposing as boundary limits {A, 7) — (0.00759, 1.35) at z = 3 and {A, 7) = 
(0.0036, 2.46) at z = 4. 

A useful way to match our adopted theoretical scheme with the observations is to 
derive from Eq. (5) the expected IGM flux attenuation (-Da) in the spectral range A G 
[1050,1170] A as a function of redshift. This can be directly compared with the observed 
deficit in quasars fluxes due to the accumulated absorption of the Lya lines in the IGM. 
By definition. 



Jl050 -f'intrl-^jaA 

Our results are summarized in Fig. ^ which collects observations from different sources. 
The data cannot be described by a single power law, since the slope of D\{z) changes in 
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Fig. 2. The deficit Da in quasars fluxes due to the accumulated absorption of the Lya 
forest. The dashed Unc is obtained by describing the IGM opacity according to the 
prescriptions by Madau (1995), while the dotted Une are those by Scott et al. (2000). 
The solid line is the model adopted in this paper. References for the data points are 
indicated. The small difference between the Madau values of Dp^ and ours at z > 4 is due 
to the contribution of Ly limit systems (Madau 1995, Eq. (13)), which was not included 
in Madau's formulation. 

the interval z e [3, 4]. The Press ct al. (1993) solution (adopted by Madau 1995) sensibly 
overestimates the IGM opacity at z < 3.5; the opposite happens with the Scott ct al. 
(2000) results when extending them beyond z > 3.5. Note, on the contrary, that our 
simple parameterization fits nicely the Da{z) evolution over the entire range of redshift, 
from z = 2 to z = 6. 

We repeated the photometric redshift estimate for the Zgpcc > 2.0 objects adopting 
our new prescriptions for the IGM absorption. For z > 2.5, only one outlier (namely 
HDF36512-1349 at Zgpoc = 3.162) is found at |zspcc - ZphotI > 0.5. In the Cohen et al. 
(2000) list the spectroscopic redshift of this object is reported as (very) uncertain. 
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Fig. 3. Comparison between spectroscopic and photometric redshifts at Zgpcc > 2. Panel 
a: the IGM opacity follows the prescription by Madau (1995); panel b: by Scott et al. 
(2000); panel c: the formulation discussed in this paper. The solid line is for Az ~ 0.0. It 
should be noted that when the Scott et al. (2000) results are adopted, the IGM opacity 
is underestimated at z > 3.5, and, as a consequence, photometric redshifts overestimate 
the spectroscopic ones. 

Excluding this source, the dispersion decreases to — 0.129 (cf. Fig. || c). More 
importantly, the tendency to underestimate galaxy redshift almost disappears (Az = 
0.042 at Zspcc > 2.5). This is not true however in the range z e [2, 2.5], where the results 
are still problematic. 

4.2. The ISM opacity and the bump at 2175 A 

In the redshift range Zgpcc & [2, 2.5], the 912 A Lyman break has not yet entered the U300 
band and therefore galaxies can still be detected in this filter, although with a major 
dimming effect by the IGM absorption. 

The U30Q magnitude does not effectively constrain z in the SED fitting method, 
because of its poor signal-to-noise ratio |^ The merit function is mainly modulated in 
this redshift range by the optical information coming from the -B450, Vgqq and /814 bands. 
It is therefore plausible that any bias in the z estimate is linked in some way with the 
^450 - VeoG and Veoe - 4 14 colors. 

At z ~ 2 the i?45o band is probing the rest-frame emission at about 1500 ± 300 A, 
while the Vgoe band samples the wavelength range 2000 ± 300 A, and the Isi4 band that 
at 2700 ±300 A. Between 1500 and 2700 A, the Leitherer et al. (1999) starburst models 
indicate that the galaxy spectral slope steepens with increasing wavelength so that, in 
general, {B450 - Vgoe) > {Veoe - hu) in the AB mag scale. 

^ Note that, paradoxically, the same spectral signature which defines so effectively the U 
drop-out sample at z > 2 causes the [/300 information to be rather useless in the SED fitting 
method. 
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Fig. 4. Observed colors of U drop-out galaxies (triangles), galaxies with spectroscopic 
redshifts in the range z e [2,2.5] (filled squares), and model colors of starburst galaxies 
in the redshift range z G [2,3.5] (see the text for details). Dashed lines show the color 
evolution of templates for increasing redshift (from bottom to top) at E{B — V) = 
0.0, 0.2, 0.4 (from left to right). Solid lines define reddening vectors aX z — 2.0, 2.5, 3.0, 3.5 
(from bottom to top), that is, the loci of templates at fixed redshift with variable amounts 
of dust. The dotted line corresponds to -B450 — Vgoe = ^606 — -^814- 

This is confirmed by Fig. ^, where the color diagram for the 100 Myr starburst model 
with Salpeter IMF and Zq is shown. Dust reddening can only strengthen this effect: 
dust absorption reddens -B450— Vgoe more than Vgoe — Is.u (see in Fig. ^ the reddening 
vectors). Moreover, at z > 2.5 the Lyman series enters the -B450 band and the -B450 — Veoe 
color is further reddened by the IGM absorption. 

In Fig. ^ we also plotted the HDFN U drop-out sample (triangles) selected according 
to the prescriptions of Meurer et al. (1999) (C/300 — S450 > 1-3, Vqoq — hu < 0.5, 
and ^450 < 26.8). A large fraction (24%) of these galaxies displays observed colors 
(i?450 — Vgoe) < (^606 — -^814) and 3 out of the 11 sources at Zgpoc G [2,2.5] are among 
these. Note that for galaxies with spectroscopic reddshifts, the color errors are < 0.02, so 
it is difficult to explain such large deviations from models in terms of photometric errors 
only. To account in a consistent way for the deviating -B450 — Vgoe and Veae ^ Isu colors, 
one could think that a dimming in the Vgoe flux is intervening in the data, at about the 
rest- frame wavelength corresponding to the 2175 A bump of the dust extinction law in 
the Milky Way (Seaton 1979) or the Large Magellanic Cloud (Fitzpatrick 1986). 
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Fig. 5. Differential dust absorption (for E{B — V) = 0.2) according to Calzetti's law 
(solid line), a superposition of Calzetti's law with a 2175 A bump of amplitude a — 0.25 
(dashed line), and the Seaton (1979) extinction curve for the Milky Way (dotted line). 

We tried to superimpose the bump C(A) on Calzetti's law, in the form of a Drude 
profile according to Fitzpatrick & Massa (1988, see also Calzetti et al. 1994): 

„2 



C(A) = a X 



^(x2-a;2)2+7a;2 ' 
where x = {^ra~^), xq = Aq"^ 



(7) 



4.6 ^.m ^ is the central wavelength of the dust 
bump, 7 = 0.95 ^im^^ is its FWHM, and = E{B - y)/(0.935 x 0.44) is the Balmer 
optical depth, as inferred from the difference in the optical depth between the nebular 
emission lines Ha and H/3 (sec Calzetti et al. 1994; Calzetti 2000). Introducing the bump, 
the differential dust attenuation curve takes the form ^(A) — [k{X) + C(A)] x E{B — F), 
where fc(A) follows Calzetti's law, and C(A) comes from Eq. (7). 

In Fig. ||, the differential dust attenuation curve is shown (for E{B ~ V) — 0.2) 
without the bump (solid line) and with a bump of intrinsic depth a = 0.25 (dashed line). 
For comparison the Milky Way extinction curve is also shown (Seaton 1979), where the 
intrinsic depth of the bump is a ^ 0.8 {a ^ 0.4 for the Large Magellanic Cloud, see for 
example Calzetti et al. 1994, Fig. 13). 

When the 2175 A bump is superimposed on the dust attenuation law, the locus 
occupied by starburst galaxies in the (^450 — Vgoe) vs. (Vgoe ~ ^814) plane (in the range 
z G [2, 2.5]) broadens to include galaxy colors with i?45o— Veoe < Vgoe — -^814 (see Fig. ^). 
Adding the 2175 A bump {a = 0.1, 0.25) enables us to better reproduce galaxy colors in 
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Fig. 6. Observed colors of spectroscopic galaxies (squares) and model colors of starburst 
galaxies in the redshift range z G [2,2.5]. Dashed lines and solid lines have the same 
meaning as in Fig. ^ but in this case reddening in colors due to dust absorption is 
obtained superimposing on Calzetti's law a 2175 A bump with amplitude a — 0.25. 
Again, the dotted line corresponds to ^450 — Vgoe ~ ^606 ^ ^8i4- As consequence of the 
rotation of the reddening vectors, model colors are able to describe galaxy colors with 
(i?450 — Veoe) < (^606 ~ -^814)- 



the range Zspcc G [2,2.5] and significantly improves the agreement between photometric 
and spectroscopic redshift estimates. While Cz = 0.206, Az = 0.167 (without the bump), 
the values decrease when the bump is allowed to — 0.094, Az — 0.060. For the three 
galaxies (see Fig. 6) which adopt a 2175 A bump value different from zero in their fit, the 

value is, on the mean, 4 times smaller than the one obtained with no bump, confirming 
the statistical significance of our improvement. 

The question of the presence of the 2175 A bump in the dust attenuation law is 
still hotly debated in the literature. Gordon et al. (1997) claimed that radiative transfer 
effects associated with the spatial distribution of dust cannot (alone) decrease the feature 
below observational detectability. Therefore, the lack of a 2175 A bump in a starburst 
galaxy must be explained as a consequence of physical and chemical properties of dust 
grain. Reversing the argument, in a starburst with LMC-like dust, geometrical effects 
alone would certainly not prevent its possible detection. 

In contrast, Granato et al. (2000) found, with their dust model for starbursts, that 
the 2175 A bump is absent due to purely geometrical effects. In the far-UV, the global 
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attenuation in starburst models is dominated by the component associated with molecular 
clouds, where stars are born, and from which they gradually escape. Clouds have such 
large optical depths that the ultraviolet light from the stars inside is completely absorbed, 
and the shape of the attenuation curve does not contain information on the UV-optical 
properties of grains. 

On the other hand, relying on the same theoretical framework, the results of Silva et 
al. (1998) lead us to conclude that virtually all the stars embedded in dark molecular 
clouds escape the regions of harder absorption in less than lO^yr (cf. also Buzzoni 2001 
on this subject). In such cases, the cirrus component dominates and a (weak) 2175 A 
bump should eventually appear in the dust attenuation law. 

Calzetti et al. (1994) did explicitly avoid the 2175 A bump wavelength region when 
calibrating their starburst attenuation law, as the possible presence of absorption bump 
would destroy a linear correlation between the UV power law slope /3 and Tg. When 
fitting a posteriori the residuals between observed spectra and the appropriate linear law 
(logF(A) = /31og A + const) according to the Fitzpatrick & Massa (1988) parameteriza- 
tion, they obtained a mean value a — 0.07 for the 2175 A bump amplitude. It should be 
noted, though, that among the objects in their sample with Tg > 0.25 {E{B — V^) > 0.1, 
see again Calzetti et al. 1994, Fig. 13) - that is, those where the presence of the feature 
is easier to detect- about 35% are compatible with a = 0.25. 

Our results support the presence of the 2175 A bump in some galaxies observed in the 
spectroscopic sample: values of a up to 0.25 allow a better description of galaxy colors 
and redshifts. Higher-quality spectra for these controversial objects are needed to assess 
beyond any reasonable doubt the actual presence of the bump. 

5. Updating our results 

The inclusion of the refinements discussed above in the SED model libraries significantly 
improved the agreement between photometric and spectroscopic redshifts at Zgpoc > 2.0. 
The dispersion decreases to = 0.117, while the systematic error is Az = 0.048 (see 
Fig. 1^), with no catastrophic outliers. 

Using the same sample in the redshift range Zgpcc > 2.0, Fontana and collaborators 
obtained = 0.232, Az = 0.144 (with one catastrophic outlier), Fernandez-Soto et 
al. (1999) obtained az = 0.299, Az — 0.186 (with two catastrophic outliers), while for 
Fernandez-Soto et al. (2001) CTz — 0.234 and Az — —0.136 (with one catastrophic outlier). 
With a smaller spectroscopic sample, Bolzonella et al. (2000) estimated ^ 0.36 (with 
one catastrophic outlier). 

A further small improvement can be obtained by taking into account dust opacity in 
the Galaxy along the line of sight to the HDFN. Introducing the effects of dust absorption 
on model fluxes in the observed frame following the Milky Way extinction curve (Seaton 
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Fig. 7. Comparison between spectroscopic redshifts and photometric estimates obtained 
with the BUZ Ubrary plus Leitherer et al. (1999) models in the entire redshift interval 
z E [0,6]. At Zgpcc > 2, the ISM and IGM opacities follow the rules introduced in this 
work. The solid line is for Az — 0.0. 



1979) with E{B - V) ^ 0.012, as inferred from the maps of Schlegel et al. (1998), we 
obtain minor changes in the redshift estimates, though in the direction of improving 
the agreement with spectroscopic redshifts. Using the BC code, for example, we find 
fJz = 0.071, = 0.002 at Zspcc < 1.5, and = 0.117, 0.036 at z^pcc > 2.0. 

Dawson et al. (2001) recently published 12 new redshifts within the HDFN. Adding 
this subset to our sample (but excluding galaxies 2-600.0 and 4-236.0, which have a very 
tentative redshift), we obtain (Jz/{1 + ^^spcc) ~ 0.040 on the largest deep spectroscopic 
sample available up to date (155 objects). This is a factor of 2.0 — 2.5 better than previous 
estimates in the literature. We find (t^ = 0.070, Az — 0.001 at Zgpcc < 1-5, and cr^ — 0.116, 
Az = 0.034 at Zspcc > 2.0. Figure ^ shows the distribution of residuals between Zgpcc and 
•Zphot, normalized to (1 -I- Zspcc)- It suggests that all estimated redshifts are confined in 
the range —0.1 < (zspcc — ^phot)/(l + ^^spoc) < 0.1, confirming the reliability of the SED 
fitting procedure. 
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Fig. 8. The distribution of residuals between 2;spec and Zphot, normalized to (1 + -Zspec)- 
5.1. The faint limit 

Our results state the accuracy of the SED fitting method when apphcd to the bright 
spectroscopic sample. When galaxies are fainter and photometric errors increase the 
accuracy obviously decreases. 

We used a simple Monte Carlo method to separate the effects of observational uncer- 
tainties from the effects due to model incompleteness ( "template mismatch" ) in defining 
the distribution of residuals between photometric and spectroscopic redshifts. We gen- 
erated 100 copies of the spectroscopic sample, adding to each galaxy flux a correction 
randomly extracted from the Gaussian distribution of its photometric error. 

We compared, for each object, the redshift obtained from the measured photometry 
Zphot and that from the bootstrapped galaxy photometry Zsim- Due to photometric errors, 
the dispersion value is ^"'/(l + Zphot) ~ 0.031, and the systematic error goes to zero. 

To quantify the contribution of model incompleteness (and only that) to the disper- 
sion, we assume that the two sources of error are uncorrelated. In this case, 



and we obtain o-^°'"p/(1 + z^pcc) ^ 0.025. 

Since the residuals introduced by template incompleteness are negligible compared to 
the photometric redshift instability caused by observational uncertainties in the interme- 
diate and faint galaxy sample (see Paper 1, Table 2), we can conclude that, at the faint 
limit, Monte Carlo methods are sufficient to estimate errors in the SED technique. 
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Fig. 9. The redshift distribution of HDFN galaxies 

In Fig. ^ we show the redshift distribution of HDFN galaxies obtained by applying 
our photometric redshift code to the entire catalog of Fernandez-Soto et al. (1999). The 
well known peak at z ~ 0.9 is clearly visible and there is a gradual decline at larger 
redshifts, with a tail at very high redshift z > 4.5. 



6. Discussion and conclusions 

We have explored the accuracy of the SED fitting method to estimate photometric red- 
shifts by comparing its results to spectroscopic measures from the deepest and most 
complete redshift catalog to date in the HDFN (Cohen et al. 2000) complemented by 
new HDF South redshift data. We used as photometric dataset the Fernandez-Soto et 
al. (1999) catalog for the HDF North, and for the HDF South the catalogue from the 



Stony Brook's group (web page http://www.ess.sunysb.edu/astro/hdfs). The libraries of 
galaxy colors were supplied by the theoretical codes of Buzzoni (2001), Bruzual & Chariot 
(1993), and Fioc & Rocca-Volmerange (1997) to reproduce normal galaxy SEDs for dif- 
ferent morphological types, and by the theoretical code by Leitherer et al. (1999) to 
describe starburst SEDs. 

A very good agreement between spectroscopic and photometric redshifts is found 
at Zgpcc < 1-5, independently of the template set used. Model libraries include enough 
information to properly reproduce galaxy SED as a function of redshift from the ultra- 



violet to near-infrared. At Zspcc > 2 the usual recipes for quantifying ISM and IGM 
absorption yield photometric redshift estimates that are systematically smaller than the 
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spectroscopic ones (A2: = 0.156), while the dispersion increases to ctz = 0.177, confirming 
previous studies. 

We focussed our attention on the prescriptions used to describe photon absorption in 
the IGM and ISM, since these processes have been shown (Paper I) to be the deciding 
factors in the redshift estimate at Zspec > 2. 

In the redshift range z € [1-5,3] the Madau (1995) prescription (adopting the Press 
et al. 1993 values) overestimates IGM opacities with respect to the values obtained by 
Scott et al. (2000) on a larger quasar sample. The observed deficit in quasar fluxes due 
to the cumulative absorption by the Lya lines of the IGM cannot be described with a 
single power law because the slope of Da{z) changes in the interval z G [3, 4]. We adopted 
the Scott et al. (2000) values at z < 3 and Press et al. (1993) values at z > 4, and a 
linear interpolation between these two estimates in the range z G [3,4]. In this way at 
z > 2.5 the tendency to underestimate galaxy redshifts disappears and at the same time 
the dispersion decreases (0-5, = 0.129). 

In the range Zgpec G [2, 2.5], 3 objects in the spectroscopic sample show optical colors 
(-B450 — Veoe < Veoa — Isu) incompatible with intrinsic model colors for starburst galaxies 
and with the dust reddening absorption law of Calzetti (1999). The common feature 
is a flux dimming in the Veoe band, which samples the restframe wavelength region 
A ^ 2000 ± 300 A at z ~ 2. Adding to the dust attenuation law a weak (a = 0.1, 0.25) 
2175 A bump, enables template colors to reproduce the observed colors of high redshift 
starburst galaxies in the {B450 — Veoe) vs. (Veoe — Isu) plane. 

When these improvements are implemented, the dispersion between photometric and 
spectroscopic redshifts at Zgpec > 2.0 becomes — 0.116, and the systematic error Az = 
0.034. In the entire redshift interval z S [0, 6], we obtain the value CTz/(1 + Zspec) ~ 0.04. 
This is an improvement by a factor ~ 2.0 — 2.5 over previous estimates in the literature. 
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